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SUMMARY 

-  Optimum  beamformers  with  a  single  look  direction  constraint 
can  suffer  from  signal  suppression  problems  when  the  optimum 
weights  are  claculated  from  the  inverse  of  the  signal-plus- 
noise  cross-spectral  matrix.  Signal  suppression  occurs 
when  the  beam  steer  direction  does  not  exactly  correspond  to 
the  signal  direction  and  this  can  occur  if  the  number  of 
fixed  beams  is  small.  The  use  of  multiple  linear  constraints 
upon  the  optimum  weights  reduces  this  signal  suppression. 
Multiple  directional  constraints  can  lead  to  ill-conditioned 
equations.  However  it  is  shown  that  the  limiting  solutions 
of  multiple  directional  constraints  are  multiple  derivative 
constraints  and  these  do  not  lead  to  ill-conditioned  equations. 
The  ability  of  various  derivative  constraints  to  prevent  signal 
suppression  is  analysed  quantitatively. 
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1.  INTRODUCTION 

Optimum  beamforming  with  multiple  linear  constraints  is  now  a  well-established 
technique  in  array  processing.  In  the  simplest  case,  a  single  constraint  is 
imposed,  viz.  unity  signal  power  response  in  the  beam  steer  direction;  the 
weight  vectors  are  then  calculated  by  minimising  the  beamformer  output  power 
subject  to  this  constraint.  This  first-generation  optimum  beamformer  is  also 
optimum  in  the  sense  that  it  (a)  maximises  the  output  signal -to-noise  ratio  (SNR) 
and  (b)  gives  a  least -mean-square  (LMS)  estimate  of  the  signal (ref . 1) .  These 

beamformers  are  usually  designed  under  the  assumptions  of  plane-wave  signals  and 
an  ideal  propagation  medium.  In  actual  operating  conditions  these  ideal 
assumptions  do  not  hold,  and  signal  suppression  can  arise  from  causes  such  as 
beam  steer  angle  errors,  phase  errors  in  the  beamformer,  and  multipath  propagation. 

When  there  are  only  a  finite  number  of  beams  to  span  the  total  bearing  angle,  any 
signal  that  is  not  exactly  aligned  with  one  of  the  beam  steer  directions  will  be 
regarded  as  an  unwanted  interference  by  the  beamformer  and  therefore  will  tend 
to  be  suppressed.  To  overcome  this  problem,  it  is  desirable  to  broaden  the 
signal  acceptance  angle  (ie  the  width  of  the  main  beam)  whilst  preserving  the 
beamformer's  ability  to  reject  interference  from  directions  outside  this 
acceptance  angle.  There  are  two  ways  of  achieving  this,  viz  .(i)  imposition  of 
a  magnitude  constraint  on  the  weight  vector  norm,  and  (ii)  imposition  of  multiple 
linear  constraints  so  as  to  preserve  the  main  beam  shape.  The  former  method  also 
improves  the  robustness  of  the  beamformer  to  phase  errors,  whereas  the  latter 
method  does  not.  This  report  is  only  concerned  with  the  prevention  of  signal 
suppression  arising  from  beam  steer  angle  errors  by  the  use  of  multiple  linear 
constraints,  and  in  particular,  investigates  the  choice  of  these  constraints. 

The  work  reported  is  part  of  a  continuing  programme  of  research  of  techniques  of 
signal  processing  for  underwater  detection  carried  out  under  Task  No.  DST  79/069. 


2.  OVERVIEW 


2,1  Notation 


Denote  by  x  =  |xj|  the  column  vector  representing  the  output  of  the  K  sensors 

at  some  frequency  f  (and  corresponding  wavelength  X) ,  The  cross-spectral 
matrix  R  of  the  sensor  outputs  is  defined  by* 


R  =  < 


x  X 


«> 


(1) 


Denote  by  z  the  K-vector  of  beamformer  weights.  The  average  power  out  of  the 
beamformer  is  then  given  by 


P  =  <z  XX  z>=  z<xx  >z  =  z  Rz 


(2) 


For  a  unit  magnitude  plane-wave  signal  incident  upon  the  array,  the  cross - 
spectral  matrix  R  is 


R  =  u  u 


(3) 


*H  denotes  the  conjugate  transpose  of  either  a  vector  or  a  matrix. 
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where  u  is  the  signal  vector  at  the  K  sensors; 

=  exp(i27r  £  t  p  ,  (4) 

T  ^  being  the  time  of  arrival  (relative  to  an  arbitrary  reference  point)  of  the 
signal  at  the  j-th  sensor.  In  general 


R  =  Q  +  u  (5) 

where  Q  is  the  normalised  noise  cross-spectral  matrix  and  is  the  signal -to- 
noise  ratio. 

2.2  Beamformers 

For  an  unshaded  conventional  beamformer,  the  weight  vector  z  takes  the 
form  ^ 


z 


c 


(6) 


where  v  is  the  so-called  steering  vector; 

V  =  exp  t  i{  j  ,<!>)] 


(7) 


where  ^^(^,0)  is  the  phase  of  a  plane-wave  signal  received  from  the  steered 

direction  (^,0)  by  the  j-th  sensor.  If  the  conventional  beamformer  is 
steered  in  the  direction  of  a  signal,  the  signal  power  p  out  of  the  beamformer 
is  ^ 


H  H  1|  H  ,2  , 

p  =  Z  U  U  Z  =  U  ul  =  1  ,  (8) 

and  so  the  conventional  beamformer  has  unity  response  in  the  look  direction. 

The  single  constraint  optimum  beamformer  uses  a  weight  vector  z  given  by 
(ref. 1,2)  ° 


Zq  =  (7  ‘  v/v*^  O' *  V  ,  (9) 

where  Q  is  the  noise  cross-spectral  matrix.  It  is  easily  shown  that  this 
optimum  beamformer  has  unity  response  in  the  look  direction  v.  A  unity  look 
direction  response  is  actually  a  single  linear  constraint  on  the  weight  vector 
z  of  the  form 


1 


(10) 


-  3  - 
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Equation  (9)  can  be  derived  by  either 

(i)  minimising  the  noise  power  z  Qz  subject  to  the  linear  constraint  (10 j, 
or 

(ii)  maximising  the  array  gain  subject  to  the  constraint  (10),  or 


(iii)  finding  the  weight  vector  that  gives  an  LMS  estimate  of  a  unit 
strength  signal  in  the  look  direction. 


All  three  methods  give  the  same  solution(ref . 1,2)  and  if  the  steering  vector  v 
and  signal  model  vector  u  are  the  same,  then  the  array  gain  g  of  the  optimum 
beamformer  is  ^ 


g 


o 


H  ^  1 
u  (7 


u  . 


(Hj 


3.  THE  SIGNAL  SUPPRESSION  PROBLEM 
3.1  Theoretical  Development 

In  a  number  of  practical  situations,  the  cross-spectral  matrix  Q  of  the  noise 
alone  cannot  be  estimated  since  the  signal  components  cannot  be  extracted  from 
the  observed  cross-spectral  matrix  R.  Thus  the  optimum  weight  vectors  must  be 
calculated  from  the  inverse  of  the  signal -plus-noise  cross-spectral  matrix. 
Furthermore,  any  plane-wave  u  not  incident  from  the  chosen  steer  direction  \ 
is  treated  as  an  interference  by  the  optimum  beamformer  and  will  tend  to  be 
rejected.  Thus  if  the  chosen  steering  direction  does  not  exactly  match  the 
signal  direction  then  signal  suppression  occurs.  This  effect  has  been 
investigated  by  Cox(ref.3)  who  derived  the  results  of  this  section. 

For  any  weight  vector  z  and  a  signal  model  vector  u,  the  array  gain  is 


g 


H 


z 


u 


H 

z 


Qz 


(12) 


Consider  first  the  ideal  case  where  the  cross-spectral  matrix  Q  of  the  noise 
alone  can  be  estimated.  For  a  given  steering  direction  (not  necessarily  the 
signal  direction)  with  associated  steering  vector  v,  the  optimum  weight 
vector  is  given  by  equation  (9),  Substitution  of  equation  (9)  into 

equation  (12)  gives  the  gain  of  the  single  constraint  optimum  beamformer  as 


®uv  I  ^ 


Su 


(13) 


where 


gjj  =  u  (f  ‘  u,  =  V  (T  ‘  V, 


(14Ca)) 
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*  H  ^  I 

=  gvu  =  ^  V 


(14(b)) 


The  notation  is  used  here  to  indicate  that  the  weight  vector  is  calculated 

from  the  inverse  of  the  noise  only  cross-spectral  matrix.  The  first  term 

in  equation  (13)  is  the  gain  when  u  and  v  are  aligned  (ie  the  beam  is  steered 
in  the  direction  of  the  signal),  and  is  the  same  as  that  given  by  equation  (11). 
The  second  term  in  equation  (13)  expresses  the  loss  in  gain  due  to  the  mismatch  he 
between  u  and  v  and  is  independent  of  the  signal  strength;  this  term  is  the 
polar  response  of  the  optimum  beamformer. 

Now  as  discussed  above,  let  the  weight  vectors  be  calculated  from  the  inverse 
of  the  signal'plus-noise  cross-spectral  matrix  R  which  is  given  by 

R  =  Q  +a^  uu^  ,  (15) 

where  Q  is  the  normalised  noise  cross-spectral  matrix  and  a ^  is  the  signal -to- 
noise  ratio.  The  weight  vector  can  be  calculated  from 

z  =  IT '  v/v^  [T  *  V  (16) 


by  the  use  of  Woodbury’s  identity; 


I  rk  -  1  ^  ^  ^  uu^Cf  * 

(Q  +  a  uu  )  =  Q  -  ^ ^ 

1  +  au  (J  u 


Substitution  of  equation  (16)  into  equation  (12)  gives  the  gain  of  the  beam- 
former  as 


^uv  I  ^ 


^u  ^v 


^uv  I  \ 


1  ^  (2a^g^ 


^u  ®v 


where  the  notation  g  is  used  to  indicate  that  the  weight  vector  is  calculated 

K 

from  the  inverse  of  the  signal-plus-noise  cross-spectral  matrix. 

If  the  steering  vector  v  is  aligned  with  the  signal  model  vector  u,  then 


^u  ”  ®v  ~  ^uv 
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and  so  g  =  g  ,  ie  there  is  no  signal  suppression  even  though  the  inverse  of 
the  signal -plus-noise  cross-spectral  matrix  was  used  to  calculate  the  weights. 

Cox  expresses  the  term  1  g^^  as  the  cosine-squared  of  the  generalised 

angle  between  u  and  v  in  a  space  where  the  inner  product  between  two  vectors 

a  and  b  is  defined  by  a  b.  Although  this  gives  a  very  elegant  method  of 

analysing  these  problems,  it  will  not  be  used  here  since  different  forms  of 
inner  product  arise  in  multiple  constraint  problems  and  hence  would  make  the 
comparison  of  the  performance  of  different  beamformers  very  difficult. 

Instead,  the  results  will  be  given  as  a  function  of  the  angular  separation 
between  two  directions. 

3.2  Examples 

The  example  used  in  this  report  is  a  horizontal  planar  array  (HPA)  and  consists 
of  25  sensors  arranged  in  5  concentric  rings,  with  a  radial  spacing  of  r 
between  adjacent  rings  and  between  the  innermost  ring  and  the  origin.  On  each 
ring  the  five  sensors  are  spaced  evenly  in  azimuth  with  one  sensor  in  the 
0  reference  azimuth  direction. 

Figures  1  and  2  illustrate  the  signal  suppression  for  beamformers  that  use  the 
inverse  of  the  noise  only  cross-spectral  matrix  and  the  signal -plus-noise 
cross-spectral  matrix  respectively  for  weight  computation.  In  both  cases  the 
ambient  noise  field  is  taken  to  be  spherically  isotropic  noise  with  uncorrelated 
noise  -30  dB  re  the  isotropic  noise.  In  figure  1,  the  second  terms  of 
equation  (13),  ie  I  g^^  ^v  plotted  as  a  function  of  the  azimuthal 

difference  between  u  and  v  (^)  at  a  frequency  corresponding  to  rA  =  0.2, 

It  can  be  seen  that  little  signal  suppression  occurs,  with  the  3  dB  loss  point 
occurring  at  ^  =  14^.  In  figure  2  the  SNR  out  of  the  beamformer  is  plotted 

as  a  function  of  ^  for  a  range  of  input  SNR  of  -15  dB  to  +6  dB.  It  is 
immediately  obvious  that  signal  suppression  becomes  more  series  as  the  input 
SNR  increases,  so  much  so  that,  for  a  given  Ada  low  input  SNR  can  lead  to  a 
larger  output  SNR  than  a  high  input  SNR.  The  aim  of  using  multiple  linear 
constraints  is  to  reduce  the  severity  of  this  signal  suppression  at  the  lower 
values  of  Ad  . 


4.  MULTIPLE  LINEAR  CONSTRAINTS 


4.1  Optimum  weights 

A  set  of  L  linear  constraints  on  the  weight  vector  z  can  be  expressed  as 

z  =  c  ,  (191 

where  V  is  a  K  x  L  matrix  whose  columns  are  the  constraint  vectors  and  c  is 
the  L-vector  of  constraint  values.  The  optimum  weights  z  are  chosen  so  as 

H  ^ 

to  minimise  the  beamformer  total  output  power  z  Rz  subject  to  the  constraints 
(19)  being  satisfied.  The  use  of  Lagrangian  multipliers  to  solve  this  problem 
is  well  known  (see,  eg  ref. 2)  and  will  not  be  repeated  here.  The  optimum 
weight  vectors  are  given  by 


WSRL-0182-TR 


6 


=  IT  ‘  v(v“  S’  *  Vf  ‘  c  , 


(20) 


and  the  beamformer  output  power  is 


p  =  (V^  IT  ’  V)”  ’  c  .  (21) 

Now  let  R  be  the  sum  of  signal  and  noise  cross-spectral  matrices  as  given  by 
equation  (15) .  Then  the  beamformer  output  power  is 

p  =  c^CV^Vq  +  uu  ^)"  *  V)"  '  c  .  (22) 

The  repeated  application  of  Woodbury*s  identity  (17)  to  the  above  equation  then 
gives  the  result: 


u  u  0^1  c^(v^q’ '  V)" '  '  u  I  ^ 

p  =  c  (V  Q^'V)"*  c  +  - 

1+0^  *  u-a^  u^Q"  *  V  (V^Q"  *  V)"  '  '  u 


(23) 


Furthermore,  if  the  i-th  constraint  is  the  signal  model  vector  u  (ie  the  i-th 
column  of  V  is  u) ,  then  it  is  not  difficult  to  show  that 


(V^q’'V)"‘  v“cr*u  =  e. 


(24) 


where  e^  is  an  elementary  vector,  ie  a  vector  with  unity  in  the  i-th  entry  and 

zeros  elsewhere.  If  the  i-th  constraint  value  is  unity,  then  substitution  of 
equation  (24)  into  equation  (23)  gives  the  result  that 

p  =  c^(v“  (J’'V)"‘c  +  (25) 

and  from  this  it  can  be  seen  that  minimisation  of  the  total  output  power  also 
minimised  the  output  noise  power.  Thus  a  beamformer  derived  from  the  latter 
minimisation  is  equivalent  to  one  derived  from  the  former  minimisation  if  there 
is  a  unity  constraint  in  the  direction  of  the  incident  signal.  If  the  unity 
constraint  direction (s)  does  (do)  not  correspond  to  the  signal  direction,  then 
signal  suppression  will  occur  in  the  beamformer  derived  from  minimisation  of 
the  total  output  power  in  a  manner  similar  to  that  discussed  in  Section  3. 

This  topic  will  be  discessed  further  in  Section  5. 

Gray(ref.2)  has  recently  shown  that,  provided  that  same  set  of  linear  constraints 
is  imposed,  one  of  which  is  a  unity  constraint  value  in  the  signal  direction, 
then  the  following  three  different  optimisation  criteria  lead  to  the  same 
solution  for  the  optimum  weight  vector.  They  are: 

(i)  minimising  the  total  output  power, 

(ii)  a  least  mean  square  fit  to  the  signal  vector,  and 

(iii)  maximisation  of  the  output  SNR. 
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Thus  all  the  results  in  this  report  are  applicable  to  beamformers  designed 
upon  any  one  of  these  three  criteria. 

4,2  Directional  and  derivative  constraints 

In  order  to  make  the  optimum  beamformer  robust  against  beam  steering  angle 
errors,  many  authors  (eg  Vural,  ref. 4)  have  suggested  the  use  of  either  multiple 
directional  constraints  or  multiple  derivative  constraints  upon  the  polar 
response  of  the  beamformer.  The  aim  of  this  is  to  attempt  to  increase  the 
angle  over  which  the  beamformer  accepts  a  directional  source  as  a  signal  and 
does  not  reject  it  as  an  interference.  This  can  be  achieved  to  some  degree 
by  maintaining  the  main  lobe  of  the  polar  response  and  not  allowing  a  null  to 
form  until  the  difference  between  the  steer  direction  and  the  signal  direction 
becomes  substantial . 

.4  two-point  unity  directional  constraint  centred  on  the  steer  direction  6  can 
be  expressed  by 


V  =  ( V  (0  -  A) :  V  (^  +  A)]  ,  c  =  (1  1)*^ 


(26) 


where  the  v(  )  are  steering  vectors. 

Similarly,  a  three-point  unity  directional  constraint  centred  on  the  signal 
direction  B  can  be  expressed  by 

V  =  [  V  (0  -  A) :  V  (0):  v  (0  +  A)]  ,  (27) 

c  =  (111)“. 


The  beamformer  performance  is  relatively  insensitive  to  the  magnitude  of  A 
for  small  values  of  A.  (For  the  HPA  used  in  this  report  very  little 
difference  in  performance  existed 'for  values  of  A  up  to  approximately  10°  for 
r/X  less  that  0.4.)  However,  if  A  is  too  large,  very  strong  signals  in 
directions  between  the  constraint  directions  may  cause  the  beam  to  split  into 
multiple  beams  with  nulls  in  the  signal  directions.  Thus  only  small  values  of 
A  are  of  interest.  However,  when  A  is  very  small,  the  matrix  V^IT  *  V  becomes 
ill-conditioned  since  the  columns  of  V  are  approaching  linear  dependence. 

This  problem  can  be  overcome  by  noting  that  the  limit  A^O  of  a  two-point 
directional  constraint  is  a  derivative  constraint  of  the  form 


V  =  ( V(0)  :  v(0)]  ,  c  =  (10)“,  (28) 

where  (  )  denotes  differentiation  with  respect  to  6  .  Similarly,  the  limit 
A-^0  of  a  three-point  directional  constraint  is  the  multiple  derivative 
constraint : 


V  =  [v(0)  ;  v(0)  :  v(0)]  ,  c  =  (1  0  0)“ 


(29) 


(The  proof  for  the  limit  of  a  three-point  directional  constraint  is  set  out 
in  Appendix  I.  The  proof  for  a  two-point  constraint  is  similar.)  Since  the 

vectors  v(5),  v(^)  and  v(0)  are  linearly  independent,  the  matrix  will 

not  be  ill-conditioned. 
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Another  form  of  multiple  directional  constraint  was  discussed  by  Vural (ref . 5) . 
In  this  case  the  constraint  values  are  those  equal  to  the  response  (magnitude 
and  phase)  of  a  conventional  beamformer,  ie  for  three  directional  constraints 
centred  on  the  signal  direction  0,  the  matrices  V  and  c  are 

V  =  [v(0  -  A)  :  v(0)  ;  v(0  +  A)  ]  . 

and 


^  V(0  -  A)”  V(0)/K  ' 

1  H 

v(0  +  A)“  v(0yK 


(30) 


The  limiting  solution  for  A  ^ 0  for  this  form  of  directional  constraint  is  the 
multiple  derivative  solution  (see  Appendix  I): 


V  =  [v(0)  :  v(0)  :  V'(0)]  . 

"1 

c  =  v(0)“  v(0)/K 

v(0)“  v(0)/K 

(31) 

• 

The  term  v  v/K  is  purely  imaginary  and  its  magnitude  increases  linearly  with 
r/X  .  (This  term  must  have  a  zero  real  component  for  the  polar  power  response 

•41 

to  have  a  maximum  in  the  direction  0,)  The  term  v  v/K  has  a  dominant  real 
component  who  magnitude  increases  as  the  square  of  r/X  . 

Thus  since  the  limiting  solutions  of  multiple  directional  constraints  are 
derivative  constraints,  and  the  fact  that  the  latter  form  of  constraint  does 
not  lead  to  ill-conditioned  equations,  only  derivative  constraints  will  be 
discussed  further. 

4.3  Discussion  and  exeimples 

An  insight  to  the  effect  of  the  various  derivative  constraints  mentioned  above 
can  be  obtained  by  examining  the  behaviour  of  these  beamformers  in  a  highly 
directional  noise  field.  For  the  following  examples  (figures  3  to  11),  the 
noise  model  used  is  spherically  isotropic  noise  plus  uncorrelated  noise  -30  dB 
re  the  isotropic  noise  plus  a  directional  source  10  dB  above  the  total  noise 
power.  The  azimuth  of  the  source  is  180  and  r/X  =  0.2.  All  the  beam- 
formers  have  one  unity  look  direction  constraint,  so  the  numerator  of  equation 
(12)  becomes  unity.  Substitution  of  equation  (20)  (with  R=Q)  into  equation 
(12)  then  gives  the  maximum  array  gain  as 

=  l/c“(vV'vr‘c  .  (32) 

Note  that  this  is  the  reciprocal  of  the  beamformer  noise  power  output  which  is 
given  by  the  first  term  of  equation  (25) . 

Figure  3  shows  array  gain  vs  asimuth  curves  for  various  constraints.  They  are 


i) 
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(i)  single  look  direction  constraint  (Type  1), 

(ii)  look  direction  constraint  plus  zero  first  derivative  constraint 
(Type  2), 

(iii)  look  direction  constraint  plus  zero  first  and  second  derivative 
constraints  (Type  3) ,  and 

(iv)  look  direction  constraint  plus  constraints  on  the  first  and 

second  derivatives  equal  to  the  conventional  beamformer  response 
(Type  4). 

The  increase  in  the  null  width  of  Types  2,  3  and  4  beamformers  over  Type  1 
beamformer  is  quite  marked,  especially  for  Types  3  and  4.  Also  of  note  is 
the  loss  in  gain  in  directions  away  from  the  source  for  the  Type  3  beamformer 
(see  explanation  below) .  The  effect  of  the  various  types  of  constraints 
can  be  easily  seen  in  the  beamformer  polar  response  plots  shown  in  Figures  4 
to  9. 

Figure  4  shows  the  polar  diagram  of  a  Type  1  beamformer  with  a  single  unity  lool 
direction  constraint  at  190  azimuth.  The  beamformer  has  a  very  deep  null  in 
the  direction  of  the  interfering  source  at  180  azimuth. 

Figure  5  shows  the  effect  of  an  extra  zero  first  derivative  constraint  (Type  2) . 
The  null  at  180  is  not  as  deep  and  a  large  sidelobe  has  appeared  at  160^. 

Figure  6  shows  the  effect  of  zero  first  and  second  derivatives  (Type  3) .  The 
polar  response  at  190  is  now  a  point  of  inflection,  and  two  large  sidelobes 
have  appeared.  The  null  depth  at  180^  is  now  only  about  4  dB.  Figure  7  shovx\^ 
the  polar  plot  for  the  same  beamformer  with  a  steer  azimuth  of  300°.  The  very 
broad  main  beam  will  lead  to  a  loss  in  array  gain  against  isotropic  noise, 
and  this  is  the  cause  of  the  gain  loss  noted  above. 

Figures  8  and  9  show  the  polar  plots  of  a  Type  4  beamformer  for  steering  azimutl. 
of  190  and  300  respectively.  The  response  of  a  Type  4  beamformer  at  azimuths 
near  interfering  sources  is  similar  to  that  of  a  Type  3,  but  at  azimuths  au\i>' 
from  interfering  sources,  the  Type  4  has  a  much  smaller  beamwidth  (at  this 
frequency  anyway)  and  so  the  loss  in  gain  against  isotropic  noise  will  be  much 
smaller. 

This  loss  in  gain  against  isotropic  noise  can  be  seen  more  clearly  in 
Figure  10  which  shows,  for  a  beam  steered  in  the  signal  direction,  the  array 
gain  in  isotropic  noise  (with  -30  dB  of  uncorrelated  noise)  vs  r/X  for  Types 
1,  2,  3  and  4  beamformers. 

There  is  no  measurable  difference  in  gain  between  Type  1  and  Type  2  beamformers 
The  lower  curve  shows  the  loss  in  gain  due  to  the  very  large  beamwidth  of  the 
Type  3  beamformer.  At  high  frequencies,  the  Type  4  gain  approaches  that  of  a 
Type  1  and  2,  whilst  at  low  frequencies  its  gain  is  similar  to  that  of  a  Type  3. 
This  is  to  be  expected  since  the  beamwidth  of  a  conventional  beamformer 
increases  as  r/X  decreases  and  hence  the  beamwidth  of  the  Type  4  beamformer 
will  increase  with  decreasing  r/X.  An  alternative  explanation  is  that  the 

constraint  values  v  v/K  and  v  v/K  in  equation  (31)  approach  zero  as  r/X  decrease 
and  hence  the  performance  of  a  Type  4  beamformer  will  approach  that  of  a  Type  3 
as  r/X  approaches  zero. 
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5.  SIGNAL  SUPPRESSION  IN  CONSTRAINED  BEAMFORMERS 
5 . 1  Theoretical  development 

The  analysis  of  signal  suppression  which  was  done  in  Section  3  can  easily  be 
extended  to  multiple  constraint  beamformers. 

For  a  noise-only  cross-spectral  matrix  Q,  the  optimiora  weights  are  given  by 
equation  (20)  with  R=Q.  Substitution  of  this  weight  vector  equation  into 
equation  (12)  gives 

I  c“(v”(T  '  Vf  '  '  u  I  ^ 

Q  c^(vV*Vr'c  (33) 

If  a  constraint  is  the  signal  model  vector  u  and  the  corresponding  constraint 
value  is  1,  then  from  equation  (24)  the  numerator  of  equation  (33)  is  unity  and 
so  equation  (33)  reduces  to  equation  (32)  as  expected. 

Let  the  signal-plus-noise  cross-spectral  matrix  R  be  given  by  equation  (15). 
Then  substitution  of  equation  (15)  and  equation  (20)  into  equation  (12)  and 
repeated  use  of  Woodbury’s  identity  gives; 

I  c“(v“(J‘V)"‘v“q’‘u 

g  ,  =  - - 

R  (1  +  o^7)^c^(v\'v)"'c  +  o'»7(c^(V^(f  (34) 


where 


7  =  -  uV  *V(vV  V)"  ‘  vV ‘u  .  (35) 

Equation  (34)  can  be  written  in  a  more  compact  form,  viz.. 


gR 


(1 


+ 


Sq 


(36) 


Again,  if  a  constraint  is  the  signal  model  vector  u  and  the  constraint  value 
is  1,  then  7=0  and  equation  (36)  reduces  to  equation  (32). 

If  the  constraint  matrix  reduces  to  the  single  constraint  of  Section  3  (ie  V=v) 
then 


«uv  1^ 


■ 


(37) 
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Also 


=  • 

Hence  from  equation  (36) , 

=  Sq 

II  *  a^Cgu  -  g^)V  ^  a^gQ(g^  -  gp) 


1  +  (2  ■  V  ’ 


which  is  just  equation  (18) . 

5.2  Discussion  and  examples 

In  the  following  examples,  the  noise  field  is  spherically  isotropic  noise  plus 
uncorrelated  noise  -30  dB  re  the  isotropic  noise  and  rA  =  0.2. 

Figure  11  illustrates  the  signal  suppression  for  Types  1,  2,  3  and  4  beamfor- 
mers  whose  weights  are  calculated  from  the  noise-only  cross-spectral  matrix. 

The  expression  in  equation  (33)  is  plotted  as  a  function  ofAO  and  has  been 
normalised,  (ie  Bq/Bq  is  plotted).  The  plots  illustrate  the  increase  in 

beamwidth  of  a  Type  3  over  a  Type  4  and  of  a  Type  4  over  Types  1  and  2.  The 

difference  in  beamwidth  between  Types  1  and  2  is  not  significant. 

For  the  signal -plus-noise  case  the  use  of  equation  (34)  provides  a  method  of 
quantitatively  comparing  the  use  of  derivative  constraints  to  prevent 
signal  suppression.  Figures  12,  13  and  14  show  plots  of  output  SNR  from 
Types  2,  3  and  4  beamformers  respectively  for  a  range  of  input  SNR  of  -IS  dB 
to  +6  dB.  These  figures  should  be  compared  with  figure  2  which  shows  the 
results  for  a  Type  1  beamformer.  It  can  be  seen  that  all  three  forms  of 
derivative  constraint  prevent  the  very  sharp  decrease  in  output  SNR  that 
occurs  in  a  Type  1  beamformer  for  large  values  of  input  SNR  and  small  values 
of  A  6 ,  The  prevention  of  sigiial  suppression  in  these  regions  is  a  highly 
desirable  property  of  any  type  of  beamformer.  For  A  6  equal  to  zero,  the 
Type  3  beamformer  shows  a  loss  in  output  SNR  compared  with  the  Type  2 
beamformer  and  this  was  noted  in  the  previous  section;  however  the  Type  3 
does  prevent  signal  suppression  for  much  larger  values  of  A  0 .  The  differences 
between  Types  3  and  4  are  negligible  except  that  the  Type  4  does  not  suffer 
the  gain  loss  of  the  Type  3  at  this  value  of  r/X  ;  for  this  reason  the  use  of  a 
a  Type  4  beamformer  is  preferable  to  the  use  of  a  Type  3, 
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6.  CONCLUSIONS 

It  has  been  shown  that  the  use  of  multiple  linear  constraints  reduces  signal 
suppression  when  using  optimally  weighted  beams.  Derivative  constraints  are 
preferable  to  directional  constraints  since  the  latter  can  lead  to  ill-conditioned 
solutions.  The  use  of  weight  vector  norm  constraints,  either  alone  or  in 
combination  with  linear  constraints,  has  yet  to  be  investigated. 
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APPENDIX  I 

LIMITING  SOLUTION  FOR  DIRECTIONAL  CONSTRAINTS 

The  steering  vectors  for  three  constraint  directions  can  be  approximated  by  the 
first  three  terms  of  a  Taylor  series,  ie 


V  (6  -  A)  =  v(0)  -  V  (6)  +  A*  V  (0 )  , 

v(0)  =  v(6)  , 

V  (0  +  A)  =  V  (0 )  +  V  (0  )  +  A*  V  (0 )  . 

These  equations  can  be  rewritten  as  the  matrix  equation; 


(I.l) 


V  =  VD 


(1.2) 


where 


D  = 


1  1  1 
-A  0  A 

A*  A* 

0 


,  V  =  I  v(0)  ;v  (0)  ;v(0)] 


(1.3) 


The  optimum  weight  vectors  are  given  by 

=  tr  ‘  V(v"ir  ‘  V)"  ‘  c 
Substitution  of  equation  (1.2)  into  equation  (1.4)  gives 

=  IT*  V^D(D^W‘VD)"‘  c 

=  IT'  VC^r'V)"'  (D*^’c)  . 


(1.4) 


(1.5) 


This  latter  equation  is  the  equation  for  the  weights  of  a  derivative  constrained 
beamformer  with  constraint  values  c  equal  to 


0 

1 

o’ 

’l  ‘ 

r  =  n«-«c  = 

1 

2A 

0 

1 

2A 

1 

5 

0 

_ J 

-2 

1 

1 

1  ^ 

.0. 

(1.6) 
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When  the  constraint  values  are  those  equal  to  the  response  of 
beamformer,  c  is  given  by 


c  = 


V(0  -  A)  v(0)/K 


v(0  +  v(0  3/K 


Substitution  of  equation  (I.l)  into  this  equation  gives 


(v^v  -  Av^v  h  ^  v^v)/K 


(v^v  +  Av^+  h  A?  v^v)/K 


and  so 


nH-i 

c  =  D  c  = 


‘H 

V  v/K 


-H 

V  v/K 


conventional 


(1.7) 


(I. 8} 


(1.9  f 
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gure  1 . 


Signal  suppression  using  weights  calculated  from  noise 
cross-spectral  matrix;  single  look  direction  constraint 
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Figure  2.  Signal  suppression  using  weights  calculated  from  signal -plus- 
noise  cross-spectral  matrix;  single  look  direction  constraint 
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Figure  3.  Array  gain  in  a  directional  noise 

field;  multiple  constraint  beamformers 
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Figure  11.  Signal  suppression  using  weights  calculated  from  noise 
cross-spectral  matrix;  multiple  constraint  beamformers 


Figure  12.  Signal  suppression  using  weights  calculated  from  signal - 
plus-noise  cross-spectral  matrix;  type  2  beamformer 
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gure  13.  Signal  suppression  using  weights  calculated  from  signa 
plus-noise  cross-spectral  matrix;  type  3  beamformer 


Figure  14.  Signal  suppression  using  weights  calculated  from  signal- 
plus-noise  cross-spectral  matrix;  type  4  beamformer 
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